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SUMMARY 


1. Intensities of Laue-Bragg scattering from 
powdered crystals of LiH and LiD were obtained at 
four temperatures, 20 # K, 77*K, 194°K, and 303 # K, 
using a recording G-M counter spectrometer and a 
special dewar specimen mount. 

2. Relative crystal structure factors, F re ^ 
and relative atomic form factors were calculated 
from the intensities. 

3. Absolute atomic form factors were obtained 
from the relative values by use of independent data. 

4. Characteristic temperatures,(5} for LIH and 
LiD were obtained from the relative form factor data. 

5 . Temperature coefficients, B t , and amplitudes 

of vibration were calculated from the( h) values. 

6 . Electron density distributions, ^ xoo , were 
obtained for LiH from the absolute atomic form factor 
data. 

7 . Radial charge distributions, U(r), were cal- 
•ulated from the P x00 . 

8. Evidence suggestive of a change in the state 
of ionization of the crystals with changes in temper- 
atur i was found. 
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AMPLITUDES OF VIBRATION AND X-RAY SCATTERING 
OF CRYSTALLINE LITHIUM HYDRIDE AND LITHIUM 
DENTERIDE: 300°K TO 20°K 

P. M. Harris C. K. Stambaugh* 

1. INTRODUCTION 

This investigation of lithium hydride and lithium deuteride vas 
begun through interest in the thermal vibrations and distribution of 
charge in the crystal. Since force constants and bonding in LiD are 
presumably the same as in LiH, experimentally observed differences in 
the scattering powers of the atoms or ions of the lattices of LiH and 
LiD for x-rays should arise solely from the mass change. Also, such 
effects in this case should be the largest observable, Bince the 
fractional change in mass is the largest possible and the fraction of 
total charge contributing to bonding, whatever the type, is also large. 

A number of investigators have made photographic studies of the 
Laue-Bragg diffraction of r A H. All are in agreement with a polar, 

NaCl-type lattice of Li + , H” ions. Ubbelohde discussed the differ¬ 
ence in cell size of LiH and LiD, attributing it to a difference in zero 
point energy. Lonsdale ^ calculated amplitudes of vibration of the 
atoms in LiH and LiD using Ubbelohde's characteristic temperatures and 

a mean atomic weight. Except for the qualitative investigation by 
(4 ) 

Bijvoet and Frederikse' the only experimental observation of the effect 
of thermal vibration on the intensities of the diffraction maxima from 
LiH vas that of Griffith Recently Ahmed has published certain 

other results which will be discussed. 

* Major, U.S.A.F., now at AIR RESEARCH LABORATORY, CAMBRIDGE, MASS: 




2. EXPERIMENTAL 


APPARATUS 

The intensity of scattering of x-rays from powdered samples of the 
hydride and deuteride at the Bragg angles was measured by means of a G-M 
counter-6pectrometer. The spectrometer was that described by Van Horn 11 , 
except that a G-M counter tube was substituted for the ionization chamber 
and a synchronous motor was attached to th- worm gear drive. The G-M 
tube was a North American Philips number 62003, with a thin mica window. 
The pulses from the tube were Integrated and transmitted as a counting 
rate by a General Radio type 1500A counting-rate meter, the output of 
which was continuously recorded by a Brown recorder. In order to have 
the time scale of the recorder chart read in degrees and minutes of 
chamber angle, the motor drive of the spectrometer was synchronized with 
the chart drive. 

For calibration purposes a Tracerlab 3C1A Autoscaler was used with 
the same G-M tube, and with the spectrometer turned by hand to fixed 
positions. 

The x-ray source was a Machlett copper target typ- A-2 tube, with 
the primal r beam filtered by nickel foil sufficiently thick to reduce 
the Cu K£ to less than one per cent of its unfiltered value. The high 
voltage, approximately thirty-five kilovolts, was supplied by an oil- 
immersed high voltage transformer and full wave rectifier unit. The 
power supply for this unit was a war surplus four hundred cycle gener¬ 
ator, belt driven by three one-horsepower electric motors. 
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To stabilize the anode current of the x-ray tube the circuit of 
Figure 1 was used. Fluctuations of plate current were detected and 
amplified by the right hand portion of the circuit, as shown. This 
aiqilified signal produced a change in the conductance of the 6 b 4 tubes, 
causing in turn a variation of the impedance of the transformer on the 
left (used as a saturable reactor), thus changing the filament current 
of the x-ray tube. Random variation of the anode current was reduced by 
this means to about three per cent. 

All data on the intensity of scattering of x-rays were obtained 
using transmission through a compacted powder-rod mounted in a modifi¬ 
cation of the dewar described by Griffith9. Figure 2 is a sectional 
view of the modified dewar. The specimen-rods, one-eighth inch in 
diameter by one-half inch in length, were mounted in brass bloqks so 
designed that the beam could completely bathe the cross section of the 
sample without striking the brass. These blockB were located in good 
thermal contact with the copper wall of the coolant reservoir. The 
x-ray beam entered and left the dewar chamber through plastic windows. 
Thermocouples were sunk in wells at the top and bottom of each brass 
block, as close as possible to the sample ends. Thin monel tubing 
supported the inner portion of the dewar in order to reduce heat leak to 
the coolant reservoir. The space between the outer wall the reser¬ 
voir and sample mount vas continuously evacuated. Insulation an d thermal 
shielding were sufficiently good that during the runs at twenty degrees 
Kelvin the loss of liquid hydrogen was about one hundred milliliters 
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per hour. The devar was supported, with its axis coinciding with that 
of the spectrometer, by a cantilever mount. 

Provision was made for mounting a sample each of lithium hydride * 
lithium deuteride, one above the other, in the dewar, to permit a direct 
comparison of the intensities of scattering of the two samples at chosen 
temperatures by a vertical displacement of the dewar. 

PREPARATION OF SAMPLES 

For this research there were available lithium hydride prepared by 
both Maywood Chemical Company and by Metal Hydrides, and lithium 
deuteride prepared by Metal Hydrides from deuteride allocated by Oak 
Ridge National Laboratory. Of the hydride preparations the former was 
the more pure and was used for all runs. 

The crystalline material was ground to pass a four hundred mesh 
sieve and contacted in a dry box under an atmosphere of dry nitrogen, 
the drying of which was accomplished by means of a train composed of 
columns of Drierite, phosphorous pentoxide, and lithium hydride, in 
that order. The third column was provided to remove any additional im¬ 
purities which might react with lithium hydride. In order to determine 
the extent of reaction of nitrogen vith lithium hydride, heavily over¬ 
exposed Do bye patterns were obtained from samples which had remained in 
the nitrogen atmosphere for the duration of a usual sample preparation. 
No evidence of lithium nitride was found; however, small amounts of 
lithium hydroxide were detected. 
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The specimen was briquetted by means of the Jig illustrated in 
Figure 3 • The tvo halves of the die in the upper left portion of the 
drawing were pressed together to form a cylindrical cavity one-eighth 
inch in diameter. This die vas clamped in the U-shaped portion of the 
vise with the cavity parallel to the arms of the U. Screws in one arm 
of the vise held the die firmly closed. The four-hundred mesh powder was 
packed into the cavity by use of an eighth-inch steel dowel-pin forced 
in by the screw in the remaining part of the vise. 

The finished specimens were measured, and then weighed in sealed 
weighing bottles containing the dry nitrogen atmosphere, all weights 
being corrected for bouyancy. They were then returned to the dry box, 
where they were placed in the dewar specimen mount, which was then seeded 
and transferred to the spectrometer and immediately evacuated. 

MEASUREMENT OF TEMPERATURE 

The thermocouples mentioned in the description of the dewar 
specimen-mount were used to determine the temperature of the samples 
during the various runs. Due to evidence of unreliability of these 
thermocouple temperatures - especially at the lower temperatures at 
which calibration was ob^-’ined by extrapolation - the shift in Bragg 6 
due to thermal contraction was used as an additional check on sample 
temperature. 

From combined thermocouple and thermal expansion data it is evident 
that the sample temperature did not differ by more than 10° from the 

[ 

[ 
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■bath temperature. As a result, for the purposes of calculation, the 
sample temperature has been taken as the normal boiling point of the 
bath used, since the intensities of Bragg scattering are not very 
temperature-sensitive. 


3- REDUCTION OF DATA 


In the experimental operation the Brown recorder was so balanced 
that a chart reading of zero corresponded to an arbitrary counting rate 
slightly less than the lowest expected background, while full scale de¬ 
flection corresponded to a counting rate slightly greater than the high¬ 
est expected reading. The balance was different for each run. Since 
the balancing circuit was such that the relation of counting rate to 
chart reading was linear, calibration could be accomplished by deter¬ 
mining counting rates at selected points with the auto-scaler. 

If N is the counting rate at any chosen counter angle, and x is 
the chart reading at this angle, the relation is: 

N = ax + p, (1) 

where a ard p are constants for a given run, but may differ between 
runs. In order to determine a and p, the data for the strongest line of 
each run were used. The area under the curve for any peak is: 


Jk fe 

= j Nd^ ■ a J xd^ + p j d^ 
pi pi m 


= a Ax + P(^2 - fa), 


( 2 ) 


where A T the area in counts per second, A x * s the are® chart units, 
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<f> is the angular setting of the counter (i.e., the Bragg angle 2©), , nd 
^2 and ^ are the upper and lower angular limits of the peak. 

If an arbitrary base line is selected, with end points Ni and Ng 
(or X;l and X2 on the chart) the area under this line is: 

A B = | {fe - fa) (N2 + Ni) 

a \ (fa " fa) ( coc 2 + ^ + O^i + P) (3) 

“ a a Bx + ^ (fa “ fa)> 

where A Bx is the area under the line in chart units. 

Then the net area between peak and base line is: 

Anet - a l ‘ % = aCA* - A Bx )* (*0 

Thus the ratio of the net area obtained from auto-scaler data to 
the net area obtained from the chart gives the scaling factor a. Sub¬ 
stitution of from scaler data, from the chart, and a as determined 
above into equation (2) gives the additive constant p. The total areas 
for both chart dita and auto-scaler data were evaluated using Simpson's 
one-third r fLe.- 

The -.■•’lues of a and p obtained in this manner compare favorably with 
those obtained by fitting scaler and chart data by the method of least 
squares, and are much more easily calculated. The agreement of scaler 
counting rotes, converted to chart units by the inverse of equation (l), 
with the corresponding values on the charts is well within experimental 
error 

( 


7 




The observed counting rates as determined from the recorder charts 
by use of equation (l) are not true counting rates because of co^ts 
missed by the finite dead time of the G-M tube. According to Kurbatov 
and Mann,^ for a random source of radiation, the correction for dead 
time takes the form: 

Nc = N/(l -rN), (5) 

where N is the observed counting rate, the dead time, and Nc is the 
corrected counting rate. 

For a pulsating source of radiation, such as the emission from an 
x-ray tube with rectified high voltage, the correction is similar, 
according to investigations by Alexander, Kummer and Klug, 1 ^ and by 
Cochran, ^ having in place of 'f an effective dead time € which is 
characteristic of the wave form of the source and the normal dead time 
of the counter. 

The value of £ was determined experimentally by the following 
method. With the x-ray source operating in its usual condition and 
the G-M tube set at a position to receive radiation of not too great 
an intensity the following four readings were taken: 

— No absorber in the beam 

Ng - One absorber in the beam 

N^ — First absorber removed and second placed 
in the beam 

N^ — Both absorbers in the beam 


■f - 
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The corrected values of these counting rates have the relations: 

NC 2 - Nci exp - mxi 
Ncg = NC]_ exp - ^ 12*2 

Nc^ = NC]_ exp - (p^ + p 2 x 2 ^ ( 6 ) 

= Nc 3 exp - njx^ 

NC 2 /NC 1 * Ncj^/Nc^ 

where pi is the linear ahsorbtion coefficient of the ’l” the absorber 
and x^ is its thickness. 

Substitution of the modified equation ( 5 ) into equation ( 6 ) gives 
for the effective dead time: 

_ N 2 N 3 - N]N 4 _ 

= n 2 n 3 (n x + N 4 ) - N x N 4 (n 2 + n 3 ) 

By this meems € (average) was found to be h.52 x 10"^ sec. 

After calibration of the charts by use of equation (l), values of 
counting rates were read from them at twenty-minute intervals of <f>. 
These were corrected for dead time by equation ( 5 ) and the experiment¬ 
ally determined € . 

In order to facilitate drawing the base line, the corrected 
counting rates were multiplied by the corresponding reciprocal Lorentz 
and polarization factors (sin 9 sin 20 / (l + cos 2 20 )) and the results 
plotted against 20 . 

Since Nc is the intensity of radiation entering the counter at a 
given angle, the net areas under the peaks on these plots are pro¬ 
portion-' 1 . to j |F | 2 , as can be seen by an examination of the expression 
for the power diffracted into the counter by a crystalline powder^; 
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in which: 


p =. p he** A . 1 + cos 2 26 
^ 0 16 rtrAo^c* 1 sin 6 sin 20 


J IF| JV, (9) 


P<j * power diffracted into the G-M counter 
P 0 * power per unit area in the primary beam 
e a electronic charge 
m * electronic mass 
c * velocity of light 

M * number of unit cells per unit volume 
= wave length of x-rays 
h 3 height of counter slit 
r = distance from sample to geiger counter 
A * area of slit limiting primary beam 
26 = angle between incident and diffracted beam 
j 3 multiplicity, i.e. number of planes of like spacing 
£ V ■ effective volume 

F = crystal structure factor, per unit cell, i.e.: 

F =§f s exp |- 2 rti (hx g ky s +Iz 3 )j , ( 9 ) 

where summation is over all atoms of the unit cell 

f g = atomic form factor for atom " s" 

h,k,£ 3 Miller indices of reflecting plane 

x,y,z 3 coordinates of atoms in fractions of cell edge 

The atomic form factor is independent of the wave length of the 

incident radiation only if the latter is sufficiently different from 

that of the absorption edge of the scattering atoms. This condition was 

fulfilled for lithium hydride and lithium deuteride. 
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From the plots as well as from the original charts, it was evident 
that in addition to lithium hydride* there was other diffracting matte* 

* Unless otherwise stated, whereever the words lithium hydride axe used 
it is to be understood that the equivalent is true for lithium deuteride. 

present. An investigation showed that all but a few of the extra peaks 
were due to the presence of lithium hydroxide. Unidentified weak lines 
were found at 67°30', 109°00', 113°00', ll8°20', and 134°20'. They are 
presumed not to arise from the lithium hydride phase. 

The effect of the impurities was to complicate the drawing of the 
base lines, since at most angles some scattering due to the LiCH or the 
unknown suatterer was superimposed on the background. However, it was 
believed that in the region between the LiH (ill) and LiH (200) and 
between the LiOH (20l) and LiH (220) the base line could properly be 
drawn tangent to the data curve, since no overlap was evident. With 
this as a point of departure, by repeated trial and error the base lines 
were drawn so that the areas under the hydroxide peaks were proportional 

p 

to the values for .1 |F | calculated for LiCH from data found in Crystal 
Structure s-*-^. 

It was discovered that each peak (lithium hydride and impurity alike) 
was distorted on the low angle side. If each peak was drawn symmetrically 
about the maximum, assuming the high angle side to be correct, then there 
remained a contribution on the low angle side which corresponded to a 
smaller peak one-tenth the height of, and one and one-half degrees away 
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from, the maximum of the peak it accompanied. These relations led to 
the conclusion that this distortion was actually a part of the peak an' 
should he included in finding the area. 

Resolution of overlapping peaks was carried out by triel and error, 
with the condition that close neighbors have the same shape. At the 
same time, correction for small peaks corresponding to the absorption 
edge of the nickel filter was made. 

Ordinates of the curve for each peak were read at fifteen minute 
intervals, and the area under the curve found by Simpson's rule. The 
net area, or corrected "integrated intensity" 

I' = kjlFl 2 (10) 

was found by subtracting the area under the base line. Results are 
shown in Tables 1 a and 1 b. 

Before correction for absorption in the sample could be made, it 
was necessary to determine the amount of lithium hydroxide present. 

This was accomplished for each run by comparing a LiOH and a LiH line 
chosen sufficiently close together that the absorption correction was 
the same for both. 

The ratio of the "integrated intensities" of the two peaks chosen 
is: 

*i (jin 2 >i . (11) 

I 2 *2 (J|fI 2 ) 2 
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TABLE la. EXPERIMENTAL INTENSITIES FOR LiH 


(Corrected for Lorentz and Polarization factors) 


h k 1 


Temp. 

■j °K 


20 

77 

194 

303 

111 

1*79 

950 

360 

883 

200 

1558 

2626 

1300 

2792 

220 

1698 

2692 

1448 

2705 

311 

1516 

2705 

1131 

2662 

222 

721 

127I* 

692 

1192 

1*00 

599 

729 

568 

759 

331 

923 

1613 

1178 

1498 

1*20 

1165 

1882 

1128 

1620 

1*22 

686 

1102 

542 

930 

121 (LiOH) 1500 

995 

1660 

1020 


TABLE lb. EXPERIMENTAL INTENSITIES FOR LiD 



(Corrected for Lorentz and Polarization factors) 




Temp, 

°K 


h k 1 






20 

77 

194 

303 

111 

521 

1582 

1498 

1766 

200 

1399 

1*397 

4410 

4801 

220 

1^55 

1*396 

4522 

4623 

311 

ll*07 

1*271 

3574 

4293 

222 

71*8 

2241 

1781 

1831 

1*00 

372 

1121 

823 

995 

331 

767 

2303 

1856 

2035 

1*20 

935 

2883 

2292 

2152 

1*22 

561 

1633 

1240 

1137 
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From equation (8) the expression for the constant for each substance 


k j ■ f ° TTTT-sr m J 2 ^ v ■ “*/• W 

lonrAm c 

Substituting in (ll) and rearranging: 



The ratio of the weights of the two substances is: 


W 1 

W 2 = d 2 M2 

where d is the weight per unit cell. 

Thus the weight fraction (xi) is : 



w 


(15) 


The two lines chosen were LiOH (121) and LiH (220). Approximate 
values for f^, fg, and lyg were used, giving: 


|F X ! = 6.47 

j = 16 

(3|f| 2 )i - 

669.5 

|f 2 | - 2.5 

J -- 12 

(j |f| 2 ) 2 = 

75-0 

The results of cal’ 

.ations are shown 

in Table 2. 




TABLE 2. WEIGHT FRACTION OF LiOH CONTAINED IN SAMPLES 


Sample 



Temp., °K 


2o 

77 

194 

363 

Li H 

.32 

.23 

•35 

.24 

Li D 

• 30 

.13 

.12 

.12 


For the mixed sample the mass absorption coefficient is: 

(d/f ) = | x LiH (u/p ) LiH + x Li0H (d/f )LiOH_} 

= x 1 ^(u/f)Li 0H “ (d/p )l±b] + (d/p )LiH 
Griffith^ found experimentally that the mass absorption coefficient for 
LiH was 0.74. Using this and the known value 12.7 for oxygen,^ the 
coefficient for LiOH can be calculated, and in turn that for the sample. 

The latter is: 

d/f * 8.73X! + .74 (17) 

With this equation, and using the measured values of p the dens' 1 ' ty 
and R (the radius) for each sample the value of (xR was found. These are 

1 Q 

shown in Table 3 - With the latter, by interpolation in the proper tables 
for © and for yR, the value of A, the absorption correction, was obtained 
for each 7>eak. 

The "intensities", -'orrected for multiplicity and absorption, can 
be expressed by: 

I "(T,hkl) = C T 2 | F (T,hkl)| 2 » ( l8 ) 

where the subscript T implies that need not be the same for all runs. 

The difference in density would cause a small variation in C^,; however 
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> 

^ it is probable that the greatest contribution to this variation would 

come from differences in the intensity of the primary beam. 

TABLE 3* QUANTITIES USED IN DETERMINING ABSORPTION 
CORRECTIONS FOR THE SAMPLES OF LiH AND LiD 


LiD LiH 


Symbol* 


Temp., 

°K 



Temp. 

°K 


20 

77 

194 

303 

20 

77 

194 

303 

*1 

0.30 

0.13 

0.12 

0.12 

0.32 

0.23 

0.35 

0.24 

w e 

.07154 

.07057 

' .07057 

.07057 

.06702 .06770 .07621 .06770 

v z 

.0913 

.0921 

.0921 

.0921 

.104 

.144 

.107 

.114 

D 

0.784 

O.766 

O.766 

O.766 

0.644 

0.594 

0.712 

0.594 


2.63 

1.44 

1*37 

1.37 

2.28 

1.63 

2.71 

1.68 

4R 

0.42 

0.23 

0.22 

0.22 

O.36 

0.26 

0.43 

0.27 


* Key to Symbols: 

Xi = Weight fraction of LiOH D * Density of sample 

W 8 « Weight of sample 4 = Linear absorption coefficient 

V 8 = Volume of sample R = Radius of sample « .159 cm 

If equation (9) is solved for the sodium chloride type structure 
of LiH, one obtains: 

Fhkl " * ( f Li ± f H )' (19) 

where the upper sign corrssponds to structure factors for which all 
indices are even, and the lower sign to those with all indices odd. 

In both cases the structure factor is real and positive, so that the 
absolute value signs in equation (l8) may be omitted. 


1 
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From equation (l8): 


^*T,hk^ = C T F (T,hkl) 

= ^C T (*Li ± 

= f' ± f' 
T,Li T,H 


( 20 ) 


f T,Li 


T,H 


( 21 ) 


If one designates by 1"^, g the value for all indices even, and by I ^ Q 
the value for all indices odd, then for a given sin 6/X one finds: 

f 1 T,e) l/Z+ 

1 1 

The necessary values of (l"p e)2 an( ^ (I"t, 0 )2 were found by 
plotting the square roots of the corrected "intensities" against 
sin 0/^ . For each set of data, two smooth curves were drawn through 

the points, and from these the desired quantities were obtained at 
selected values of sin e/l . 

The Debye-Waller expression for the temperature-dependent atomic 
form factor^ is: 


f* « f 0 e 


-Mp 


f R 6 


„ sin 2 © 
-ST-^ 


( 22 ) 


where f , the form factor of the atom at rest, and Brp, a temperature- 
R A 

dependent quantity, are different for each type of atom. The quantity 
Bp will be discussed in greater detail later. 
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From equations (20) and (22): 



and similarly for the hydrogen. If the values are taken so that 
sin ©1 = sin.©2> then equation (23) simplifies to: 



e ^P " (B303 

Cj 



(23') 


and, taking the logarithm: 


In 




Thus an extrapolation of equation (2k), which is linear in 


(sin 0/A ) , to sin e/A equal to zero will give the value of ln - ^ 


c 303 


which is needed in determining the absolute atomic form factors. (See 
Figure 4). 

For the purpose of placing the scattering factors on an absolute 
basis, i.e., to evaluate C>r of equation 18, a direct comparison of LiH 
with NaCI by substitution was made at room temperature.* 

A plate of LiH cleaved parallel to a cube face was used for the 
measurement of the integrated reflection, ^ , of the (200) plane of LiH. 
In the case of NaCI, a large rock salt crystal was used, the integrated 


* We are indebted to Mr. Paul D. Splitstone for carrying out this 
measurement. 


18 



reflection being measured by reflection from the cube face for the (200), 
(400), and (600) orders. 

Thus: 


pLiH 

fNaCT 


„2t-« 2 1 + cos^ 20 cos 0 
N F --sTn-^ - e 


t 

cos 0 


LiH (200) 


N % 2 


1 + cos 2 20 
sin20 


NaCl (hOO) 


The values of James and Firth ( 20 ) were used for F NaC1 (uncorrected for 
extinction since no correction was made with these data). Extinction 
should have little effect, in any case, on the (600) reflection as shown 
by James and Firth's data . The values used are: 


F NaCl 
T— 


- 17.16 (200), 11.23 (400) and 6.84 (600) 



1 . 20 ; 


'LiH (200) 
’NaCl (400) 


8.49; 


1H (200) 
aCl (600) * 


21.8 


thickness (t) of LiH plate * 0.144 cm 

N 5 no. of unit cells per unit volume =_A_ 

a 0 3 

Using the NaCl (200),(400), and (600) reflections data, F^g is > 

respectively, 5-28, 5*72, 5*40 with an over-all average of 5*48 at room 
temperature). 

Using the value of 5*48 as the scattering power per unit cell for 
the (200) reflection of LiH at 303°K, the value of C303 was determined. 
Values of C^ for the individual ions at 77° and at 20°K were then ob¬ 
tained from these plots. 

As an alternative procedure, using an average value of the charac¬ 
teristic temperature, (5), of 851°( Table 5) values of f R , the ’test" 
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value of the scattering factor were calculated from the values of f^Q^ 
for both Li and H. (The values of f Li and fg were independently re¬ 
determined for this purpose by one of us from the experimental (FLiH)303 o 
values.) Then from these values of f R , values of fand fgQ 0 were 
calculated. 


Putting 


(h 2 + k 2 + l 2 ) 

- = f, 


-BA (h 2 + k 2 + l 2 ) 


then 

= £ , where fm is the relative scattering 
m A 

factor (experimental) for the ion and frp is the value as calculated 

above; m runs over nine values (h, k, ?.). 

2 

The apparent squared error € m was minimized by the least squares 
method to yield values of C T for both Li and H at 77° and 20° K. The 
values of C T obtained in this way are compared in Table 4a with the 
values read from the extrapolation plots. It is clear that, within ex¬ 
perimental error, the values are not different. 

Although the data for LiD were related through experiment with the 
data for LiH, it is clear that these also have been correctly treated 

since the values of f D for Li in LiD are prar tically identical with 
n 

those in LiH. 
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TABLE 4a. PROPORTIONALITY CONSTANTS FOR CONVERTING 
RELATIVE FORM FACTORS TO ABSOLUTE FORM 
FACTORS: LiH, (4 Ct ) -1 


Temp., 

Least 

squares fit 

Graphical extrapolation 

*K. 

Li 

H 

Li 

H 

20 

0.0641 

0.0740 

0.0639 

0.0704 

77 

0.0522 

O.O65O 

0.0524 

0.0639 


TABLE 4b. PROPORTIONALITY CONSTANTS FOR CONVERTING 


RELATIVE FORM FACTORS TO ABSOLUTE FORM 
FACTORS: (4 C-p) -1 


Temp., 
*K. 

LiH 


LiD 

Li 

H 

Li 

D 

20 

0.0639 

0.0704 

0.0627 

0.1060 

77 

0.0524 

0.0639 

0.0408 

0.0680 

303 

0.0503 

0.0503 

0.0378 

0.0566 


The exponent B^ in the Debye-Waller expression can be written 2 - 1 -: 


u = 6h 2 

mk <5)i 


4 (x) 


1 

+ IT. 


(25) 


where 

h = Planck's constant 
m - mass of atom at lattice point 
k - Boltzman's constant 
(§) as Debye characteristic temperature 

_ hVft 

v m = maximum frequency of elastic vibration of solid 
x =w where T is the absolute temperature, °K. 


^(x) 

x 



(26) 


( 
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The function tf(x) was evaluated numerically by Debye^, and is tabulated 
in his paper for x up to 20. For larger x, Debye's equation: 


£L*1 l/«! , e ' X (x t l)-e' 2X (27) 

X 6 4 ■> 

was used. 

Using equation ( 25 ), the temperature dependent form factor may be 
written: 



( 22 ') 


so that 


l n £t_ * 6h2 f/V00) _ U(x) \ 

f303 i/\ x / 303 V x A 


TJ 




( 28 ) 


By multiplying both sides by (/£ /sin Q) 13 one obtains: 


iifee ln f^03 " 3®|^o3 - (29) 

where e)T is a constant for a given temperature and atom type. A plot 
of the functions: 


g(D 




hCD) = m k S t 
^ “5h2— 




against (u) gives an intersection at the value of (e) corresponding to 
the solution of equation ( 29 ). 

By ohis method four values of (S) for LiH and four for LiD were 
found (See Figure 6). These are presented in Table 5* 
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TABLE 4c. ABSOLUTE ATOMIC SCATTERING FACTORS 





FOR THE 

ATOMS OF 

LiH AND 

LID 



Temp., 








®K. 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 





f T (Li+) (LiH) 




20 

l7023 

0.928 

O. 85 I 

0.780 

0.711 

0.652 

0.544 

0.422 

77 

1.021 

O .922 

O .838 

0.780 

0.713 

0.644 

0.534 

0.4l4 

303 

1.003 

O .895 

0.810 

0.739 

0.664 

0.589 

0.483 

0.371 





f T (H-) 





20 

0.39 

0.26 

0.168 

0.099 

0.077 

0.071 

0.04- 

0.04“ 

77 

0.39 

0.26 

0.160 

0.102 

0.083 

0.065 

0.04- 

0.04- 

303 

0.32 

0.20 

0.122 

0.071 

0.051 

0.04- 

0.02- 

0.01 


f T (Li + ) (LID) 

20 

1.017 

O .916 

0.828 

0.765 

0.690 

0.615 

O .508 

0.396 

77 

1.007 

O .905 

0.824 

0.755 

0.682 

0.607 

0.501 

0.388 

303 

0.966 

0.853 

0.751 

O .676 

0.597 

0.510 

0.408 

0.302 






f T (D-) 




20 

0.53 

0.39 

0.22 

0.148 

0.116 

0.102 

0 . 06 l 

0 . 06 l 

77 

0.53 

0-37 

0.22 

0.142 

0.110 

0.095 

0 . 06 Q 

0.041 

303 

0.45 

O .30 

0.171 

0.095 

O.O 67 

0.061 

0.020 

0.016 
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TABLE 5- CHARACTERISTIC TEMPERATURES 
FOR LiH AND LID 


Temp., # K 

Atom type 

LiH 

T^D 

20 

Li 

850 

638 


H 

835 

636 

77 

Li 

863 

637 


H 

855 

642 

Average 


B 5 T 

U 38 

Ubbelohde^ 


815 

6 ll 


From 

the average 

characteristic 

temperatures, the temperature 

coefficients Bj of Table (6) were calculated 

; and from these and the 

absolute form factors 

of Table (4), the form factors for 

the atoms at 

rest were 

determined. 

(See Table ( 7 ) 

and figures (5a) and (5^)• 


TABLE 6. TEMPERATURE COEFFICIENTS, 




Bt IN (kXr FOR 

Li, H AND D 






Temp., °K 


Crystal 

Atom type 

0° 

20° 

' 77 b 

303* 

LiH 

Li 

0.485 

0.487 

0.510 

O .836 


H 

3-35 

3.36 

3.52 

5.78 

LiD 

Li 

0.648 

0.653 

0.710 

1.38 


D 

2.21 

2.22 

2.42 

4.71 


TABLE 7- AVERAGE FORM FACTORS, f R , FOR 






Li, H, 

AND D ATOMS AT 

rw * 

REST 








Sin e/X 




Atom type 

0.25 

0.30 

0.35 

0.40 

0.45 

0.50 

0.55 

0.60 

LiH: 

Li 4- 

1.055 

•970 

.901 

.849 

.788 

.731 

.629 

.499 

LiD: 

Li + 

1.055 

•970 

.897 

.849 

.792 

.723 

.617 

.499 

LiH: 

H~ 

.467 

.369 

.244 

.175 

.162 

.158 

.122 

.138 

LiD: 

H* 

.609 

.467 

• 304 

.203 

.179 

.179 

.122 

.102 

Mean 

H" 

• 539 

.418 

.274 

.189 

.171 

.169 

.122 

.120 
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The exponent Bj can also be written as^: 

B t « 8 it 2 u| ( 30 ) 

where Us is the component of the displacement of the atom in the direction 
of the normal to the 'reflecting' plane. For an isotropic crystal the 
total mean square displacement is related to U§ by the relation: 


U 2 = 3 U§ ( 31 ) 

Hence the total root mean square displacement or amplitude of vibration 
of the atoms may be expressed as: 

(32 > 

The root mean square displacements calculated by Lonsdale^ were 
average values expressed by: 

(33) 

2 

where M^ and Up are the mass and mean square displacement, respectively, 
of the ith kind of atom. 

The + otal root mean square displacements for the individual atoms 
sure shown in Table 8, wV.e Table 9 gives the average values for com¬ 
parison with Lonsdale's calculated quantities. 


TABLE 8. TOTAL ROOT MEAN SQUARE DISPLACEMENTS 
OF ATOMS OF LiH AND LiD 


Crystal 

Atom 


Temp., 

°K 


0 4 


77" 

W5*~ 

LiH 

Li 

.13b 

.136 

• 139 

7TW 


H 

• 356 

.358 

.366 

.469 

LiD 

Li 

• 157 

•157 

.164 

.229 


D 

.290 

.290 

• 303 

.423 
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TABLE 9- AVERAGE ROOT ME AH SQUARE 
i AMPLITUDES FOR LiH AMD LID 





Temp., *K 



Crystal 

0 

20 

77 

303 

293* 

LIH 

.179 

.100 

.184 

.236 

.24 

LID 

•195 

.195 

.203 

10 

-4 

o^ 

.29 


* Lonsdale's calculated values.8 

Fro* the atonic form factors for the llthiun and hydrogen at 20*K 
(Table 4b) the crystal structure factors for lithium hydride were de¬ 
termined by use of equation (19 )• With these the electron density p(xoo) 

pit 

along the cell edge was calculated, using the expression*^: 

oO 

Vp(xOO) F hkl c08 2jt (33) 

—CO 

Since the resulting electron density shows diffraction effects because 

of termination of the series,(Figure 7), it appeared desirable to ex- 

sis. 0 

trapolate the atomic form factor curves to larger values of —. 

At the same time it was decided that the electron density of the atoms at 
rest would be more informative. Consequently the atomic form factor values 
for lithium and hydrogen in Table 7 v«re extrapolated according to the 
method of Bouman 2 ^ . The extrapolation functions were: 

for Li fe = 1.920 exp 3*75 (sin e/A ) 2 

o (3*0 

for 1 fe - 0.455 exp 4.17 (sin e/a ) 2 
with the resulting values of F ^vi the electron density along the cell 
edge (Figure 7) was obtained. 


t 
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The lowest point of the curve in Figure 7 was chosen as the boundary 
between the lithium and the hydrogen ions. The distance of this ooint 
from the origin (o-70 A) is the same as that given by Wyckoff 2 ^ for the 
radius of the Li ion. With the calculated electron density and the 
usual expression for the radial charge distribution: 

U (r) = 4 nr 2 (r) (35) 

the curves of Figure 8 were obtained. 

The maximum in the curve for the lithium corresponding to the radius 

O 

of the K shell, occurs at 0.31 A., while a minimum in the curve for the 

O 

hydrogen occurs at about 0.5 A. 

The area under the radial distribution curve (i.e., the charge per 
ion) for the lithium is 1.49 electrons, while that for the hydrogen is 
2.00 electrons. 


4. DISCUSSION OF ERRORS 

The relative error in the experimentally observed counting rates, 
varies from about one-half per cent for the largest ccsuiting rate to 
about five per cent for the lowest, with a mean value of approximately 
two per cent. 

The correction for missed counts (equation 5) involves the use of 
an experimentally determined effective dead time which may be in error 
by as much as fifty per cent. However, even with this large an error, 
for the extreme case of the highest observed counting rate (300 cps) the 
relative error in the corrected counting rate due to the error in € is 
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six per cent, while for the lowest it is one-half per cent. Thus the 
mean relative error arising from this cause in the areas under th^ peaks 
of the corrected curves is around two per cent. 

A far more serious source of error in the 'integrated intensities” 

(I*) is the uncertainty involved in drawing the base line and resolving 
the peaks, as a result of the presence of impurity. 

Two separate spectroscopic analyses of the materials used indicated the 
presence of not more than trace quantities of such impurities as sodium, 
calcium, aluminum, and other heavy metals. 

The method used depends on crystal structure factors calculated for 
lithium hydroxide, which are based on not very accurately estimated 
parameters for the oxygen atom. Further the resolution of overlapping 
peaks is somewhat arbitrary. While no calculation of the resulting error 
is possible, it is estimated that the error in intensity estimation could 
amount to as much as five per cent. 

The relative intensities (I*) which are obtained by correcting I'for 
multiplicity and absorption are within a given run as accurate as the 
integrated intensities from which they are derived. This is true because 
although tne absorption corrections may be in error they are essentially 
independent of the scattering angle for such light absorbers. However, 
as soon as an attempt is made to place all runs on a common basis, the 
magnitude of the absorption correction becomes important. In calculating 
the latter the factors subject to error are the density and the mass 
absorption coefficient of the sample. The relative error of the volume 
of the sAmple is six per cent and that of the mass is about one-tenth 
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per cent, so that the relative error of the density is six per cent. 

The correctness of the mass absorption coefficeimt depends in on 

the accuracy of the composition of the samples, and since the latter 
is determined from the ratio of tvo lines on the I' plot it may be off 
by as much as ten per cent. Further, the experimental mass absorption 
coefficient for LiH as obtained by Griffith^ was used in the calculation, 
and since his bricquetted powder sample probably contained lithium hy¬ 
droxide (for which no allowance was made) his value is undoubtedly high 
by an unknown amount. The result of these factors is that the absorption 
correction may be in error by as much as ten per cent. 

The calculation of the characteristic temperatures, the exponents 
Brp and the amplitudes of vibration depend on the slope of equation (24) 
and on the temperature at which the data were taken. Since the slope of 
equation (24) is independent of the magnitude of the proportionality 
constant, the uncertainty of the absorption correction does not affect 
the above quantities. The function, G(@ ), used in determining (H) is 
temperature-dependent, and a five-degree temperature error would cause 
an error of four degrees in (g), an amount which is negligible compared 
to the two per cent" error that could result from a five per cent error 
in the integrated intensities. Since at the temperatures for which 
measurement, is uncertain p^x)/x is small compared to l/4 (equation 25), 
the exponent Bj is relatively insensitive to temperature variation. A 
five degree error in temperature causes a negligible error in Bp, and 
hence in the amplitudes of vibration. Thus the primary source of error 
in Bt and is the two per cent error possible in (5). 


29 





The error in the absolute form factors as tabulated is difficult to 


evaluate. The principal error should arise in the determination of the 
proportionality constants, Crp, and the fixing of the absolute value of 
one of these. 

It is believed that the absolute value of C 303 for LiH is good to 
± three per cent. Further, the individual ion-values (relative) of Crp 
appear to be good to one per cent for the case of lithium and five 
per cent for the case of hydrogen. 

Ahmed^® has recently determined the values of (FhiH)rel room 
temperature from a single crystal of LiH using a photographic method. 

His results at this temperature agree rathei iell with the corresponding 
values reported here except for the reflections ( 220 ) and (224) which 
appear to be too large by comparison (ca. fifteen per cent). 


5. DISCUSSION OF RESULTS 

It is convenient to discuss first those results which depend only 
on relative intensities and are independent of the method employed for 
placing the measurements on an absolute basis. 

A significant feature of these measurements is that the ln(f’ 303 /f * «p) 
plots for lithium and hV~ogen (Figure 4) do not extrapolate to the 
same point. According to Bquation ( 18 ) the constant Crp for each run 
should be the same for the lithium as for the hydrogen. An error in 
drawing the base line for diffracted intensities could give rise to the 
crossing point noted; however, to eliminate the error would require that 
the bes^. 1 /ine be moved up, for which there is no justification in the 
appearance of the data, and in addition this would introduce curvature 
into the plot of ln^-^/f' <p) • 
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The most plausible explanation appears in a reexamination of 
equation ( 23 ')* In "this it is assumed that the form factor for the *tom 
at rest (fg) for the lattice at room tem p erature is the same as that for 
the lattice at temperature T. If this is not true, then equation (23) 
should be written as: 



c 3Q3 (f R ) 303 

C T (fft> T 


exp - £( 1*303 - Bcj)(sin Q/X ) 1 (23'') 


(A plot of ln(f'203 /^’t) would not necessarily be linear in (sin 6/X ) 2 .) 

Nov a change in f^ between the tvo temperatures implies a change 
in the state of ionization of the atom concerned, so that if represents 
the effective charge of the ion at temperature T, then the logarithm of 
Equation (23") would extrapolate to: 


In £303i303 slnce t . as 0 . 

C T Zrj, A 

Since the total charge of the cell must remain constant, an increase 
in Zfji for the lithium with decreasing temperature would have to occur 
at the expense of a corresponding decrease in for the hydrogen. Such 
a redistribution of charge would make the constant for the lithium 
smaller than that for the hydrogen, in agreement vith the observed results. 

Added confirmation r !" this conclusion appears in the fact that if 
the ratio f' 77 /f '20 is used for the extrapolation the difference between 
the lithium and hydrogen constants becomes much smaller, which would 
imply that the change occurs between 7T°K and 303°K. 
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One can make an estimate of the sort of charge redistribution this 

(£ 3 £ 3 \ 

corresponds to: The average value of the ratick£T_/_ Li .. "or both 

/£ 30 |\ 

\ C T / H 


LiH and LiD at both (lower) temperatures. As was discussed before, this 
means that f 2 o° for H multiplied by gives a value too 


small at room temperature and vice versa for Li 

sin 0 


Now, since 


v6in9\ 

X* ) 


Z as 


A 


o, 


(Z 303 ) 

v Zr J h = l.lL. 

(230f) 

V Z T J Li 


One can proceed further if it is assumed that LiH (or LiD) is completely 
ionized at room temperature; that is, Z Li * Zjj = 2 electrons. Then at 
T (20° or 77°), Z Li = 2 + <fz, Z H = 2 - /Z. S« §- ■ — - l.lU; f Z = 0.13- 
Thus, on these assumptions, at a sufficiently low temperature LiH 
and LiD are not completely ionic, but about 87 # of the charge is trans¬ 
ferred. It is not surprising that the charge transfer is about the same 
at 77° as at 20°K, since the values of are so nearly the same and 
are so near their zero-point values. 

The possibility of appreciable covalent bonding in LiH has been 
suggested by Ewing and Seitz^. They have estimated thet- for the lattice 
at rest the bonding should be apprnxinrateiy^ fty^Yinmmflro f 35 $ coulombic. 
Ahmed^ has interpreted his room-temperature data on LiH by extrapolation 

of f Li> f H = 

f Li ^ Z Li> f H * Z H as ~ 'j - 9 — 

Using this extrapolation he obtains a value of Z L i of 2.75 - 0.25 

electrons, (and, of course, 1.25 % 0.25 for the hydrogen). Since the 
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plots of f vs sin B/X have a large and increasing (negative) slope in 


the region of small —, such an extrapolation is open to considerable 
question. 

One possible explanation of the change in the state of ionization 

could be found in the symmetry of the lattice. For the case of the 

atoms at rest, each lithium atom is coordinated with six hydrogen atoms 

on orthogonal axes, and conversely. This condition is favorable to a 

resonance of covalent bonds throughout the six bonding directions. As the 

farther 

temperature is increased, the nuclei begin to vibrate^ from their equi¬ 
librium rest position, destroying the symmetry that favored resonance 
and thus favoring a shift to coulombic binding with its greater flexi¬ 
bility. Such a process would be compatible with the suggested increase 
in ionization of the atoms with increasing temperature. 

There is theoretical justification for a change in the characteristic 
temperature 2 ^ with the temperature of the crystal; however, its only 
effect on the logarithmic plots would be to introduce a curvature which 
was not observed. Further, the characteristic temperatures shown in 

7 

Table 5 are in excellent agreement with those calculated by Ubbelohde 1 
(6l5°K for LiH and 6ll # K for LiD), especially when account is taken of 
the fact that his values are based on specific heats, and according to 
James 2 ® the calculated from x-ray measurements should be between 
three and five per cent higher because of a difference in the method of 
obtaining the mean of the maximum of frequency for the longitudinal and 
transverse waves. From these values of (h) and the relationship: 

hvn»x = *(H) (40) 


33 



one finds that the mean maximum frequencies of vibration are 1.77 x 10^ 
sec-1 f or the LiH and 1.33 x 1013 sec-1 f or the LiD. These con >spond 
to reststrahlen wave lengths of 17 ^ and 22 ^ for the LiH and LiD,, 
respectively. Thus it should be feasible to obtain a confirmation of the 
(§) values by infrared measurements. 

The amplitudes of vibrationV agree well with Lonsdale's ® 
calculations (Table 9)> & not very surprising result since she used 
Ubbelohde's characteristic temperatures. The amplitudes for lithium and 
for hydrogen are in the proper ratio for their relative masses, and the 
change in amplitude resulting from the substitution of deuterium for 
hydrogen is as would be expected. 

The interpretation of (h) } Bp and depend on the degree to which 
the crystal of LiH approximates the ideal model of the Debye theory. The 
most serious departure from this model is that, instead of a crystal of 
atoms of one kind, lithium hydride is a binary crystal. Hence, there may 
be a wide separation between the acoustic and optical vibration bands. 

From the charge density determinations one obtains for the radius 
of the lithium ion at rest the value 0.70 angstroms, and by subtraction 
from the half cell-edge, the value 1-33 angstroms for the hydrogen ion 
at rest. By apportioning the change in half cell-edge with increasing 
temperature between the two ions in the ratio their amplitudes of 
vibration, the ionic radius for lithium at room temperature remains 0-70 
angstroms, but that for hydrogen becomes 1-3^ angstroms. According to 
equation (22) these radii are for the atoms in the state of ionization 
found at room temperature. The lithium radius is in exact agreement 
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with that reported by WyCkoff" for the Li t ion, v, al'.e the radias for 
the hydrogen is within six per cent of this value i-~r the H“it.. 

The radius of the K shell for the lithium ion, as determined from 
the maximum of the radial charge distribution for lithium ie 0,31 
angstroms, which is in acceptable agreement with theoretic, calculations 
by Wilson 2 ^. (His value is 0.28 angstroms.) The total cL-jxge for the 
lithium ion (1.49 electrons) which is determined from the integration of 
the radial charge distribution ever (rj appears to be much too lew, since 
a partially ionized lithium atom should have a charge of greater than 
two electrons. 

Thus, 3*5 of the lour electrons are found in the spheres drawn 
about the Li and H centers. Since integration of the charge density 
over one-fourth of the unit, cell must yield just fn „r e_-c» T, ens as a 
result of the F 00O term in the summation, the unaccounted.-fcr charge 
is in the interstitial space r.c 4 covered fy the sph-rn a^ integration 
limits chosen. In any case, this integration is Mit j*.-1 to consider able 
error at large values of the radius. 

Considerable confidence in. the at sc.lux-? mag: it he 4 o* tie scattering 
powers is contributed by the fact that the charge : y ^ co falls 
•ic nearly to zero between the 11 r.htut. aiid the hydroger - 

Ewing and Seitz30 calculated the charge ais- ributioc. for the LiH 
lattice at rest using an approximate, self■■■-onsisiend field tsethod. 

They list no quantitative results, sc if was not possible to calculate 
scattering factors for ccmparati'■*<? purposes. Their T ! ,r) for Li and H 
is compared with the results obtained here in Figure..- B.-x, t. For the 





purposes of comparison, since the ordinate scale va:> unknovn, the scale 
has been chosen so as to give as good agreement \ih: the U(r) Mt-hium 
as possible. 
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